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Abstract 

O ' 

^-ivVe consider the double pulsar PSR J0737-3039A/B binary system as a laboratory to locally test the orbital effects induced by 
I an uniform cosmological constant A in the framework of the known general relativistic laws of gravity, and the DGP braneworld 
^jynodel of gravity independently of the cosmological acceleration itself for which they were introduced. We, first, construct the 
i_iratio R = Aui/AP of the discrepancies between the phenomenologically determined periastron rate Co and orbital period Pb 
and their predicted values from the 1PN wipn approximation and the third Kepler law P^ ' . Then, we compare its value \R\ = 
(0.3 ± 4) x 10 11 s 2 , compatible with zero within the errors, to the ratios Ra and -Rdgp of the effects induced on the apsidal rate 
and the orbital period by A and the DGP gravity; we find them neatly incompatible with R being Ra = (3.4 ± 0.3) x 10~ 8 s -2 



> 

in 



and -Rdgp = (1-4 ± 0.1) x 10 s , respectively. Such a result, which for the case of A is valid also for any other Hooke-like exotic 



force proportional to r, is in agreement with other negative local tests recently performed in the Solar System with the ratios of 
f—^ the non-Newtonian/Einsteinian perihelion precessions for several pairs of planets. 
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1. Introduction 

■ Since, at present, the only reason why the cosmological 
constant 1 A is believed to be nonzero relies upon the ob- 
served acceleration of the universe (Riess et al., 1998; Perl- 
mutter et al., 1999), i.e. just the phenomenon for which A 
was introduced (again) , it is important to find independent 
observational tests of the existence of such an exotic com- 
ponent of the spacetimc. 

In this paper we wish to put on the test the hypothesis 
that A / 0, where A is the uniform cosmological constant 
of the Schwarzschild-de Sitter (Stuchlk, 1999) (or Kottler 
(Kottler, 1918)) spacetime, by suitably using the latest de- 
terminations of the parameters (see Table 1) of the double 
pulsar PSR J0737-3039A/B system (Burgay et al., 2003). 
The approach followed here consists in deriving analyti- 
cal expressions Oa for the effects induced by A on some 
quantities for which empirical values O meas determined 
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from fitting the timing data exist. By taking into account 
the known classical and relativistic effects Oknown affecting 
such quantities, the discrepancy AO = O meas — Oknown is 
constructed and attributed to the action of A, which was 
not modelled in the pulsar data processing. Having some 
AO and Oa at hand, a suitable combination C, valid just 
for the case A ^ 0, is constructed out of them in order 
to compare C meas to Ca- if the hypothesis A ^ is cor- 
rect, they must be equal within the errors. Here we will use 
the anomalistic period Pb and the periastron precession Co 
for which purely phcnomcnological determinations exist in 
such a way that our C is the ratio of Acj to AP D ; as we will 
see, this observable is independent of A but, at the same 
time, it retains a functional dependence on the system's pa- 
rameters peculiar to the A— induced force and of any other 
Hooke-like forces. 

The present work complements (lorio, 2008) in which a 
similar test was conducted in the Solar System by means of 
the latest determinations of the secular precessions of the 
longitudes of the perihelia of several planets. The result of 
(lorio, 2008) was negative for the Schwarzschild-de Sitter 
spacetime with uniform A; as we will see, the same conclu- 
sion will be traced out of this paper in Section 2.1. 



A complementary approach to explain the cosmic accel- 
eration without resorting to dark energy was followed by 
Dvali, Gabadadze and Porrati (DGP) in their braneworld 
modified model of gravity (Dvali et al., 2000). Among other 
things, it predicts effects which could be tested on a lo- 
cal, astronomical scale. In (Iorio, 2008) a negative test in 
the Solar System was reported; as we will see in Section 3, 
PSR J0737-3039A/B confirms such a negative outcome at 
a much more stringent level. 

The conclusions are in Section 4. 

2. The effect of A on the periastron and the orbital 
period 

The Schwarzschild-de Sitter metric induces an extra- 
acceleration 2 (Rindler, 2001) 



A A = ^Ac 2 r, 



(1) 



where c is the speed of light; eq. (1), in view of the extreme 
smallness of the assumed nonzero value cosmological con- 
stant (A f» 10 -52 m -2 ), can be treated pcrturbatively with 
the standard techniques of celestial mechanics. In (Kerr et 
al., 2003) the secular precession of the pericentre of a test 
particle around a central body of mass 9JI was found to be 
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2n(°) 



where 
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GM 



(2) 



(3) 



is the Keplerian mean motion; a and e are the scmimajor 
axis and the eccentrity, respectively, of the test particle's 
orbit. Concerning a binary system, in (Jetzer and Sereno, 
2006) it was shown that the equations for the relative mo- 
tion are those of a test particle in a Schwarzschild-de Sitter 
space-time with a source mass equal to the total mass of the 
two-body system, i.e. 9Jt = triA + ttib- Thus, eq. (2) is valid 
in our case; a is the semi-major axis of the relative orbit. 

Following the approach by Jetzer and Sereno (2006), we 
will now compute Pa, i-e. the contribution of A to the or- 
bital period. One of the six Keplerian orbital elements in 
terms of which it is possible to parameterize the orbital mo- 
tion in a binary system is the mean anomaly M. defined as 
M. = n(t — Tq), where n is the mean motion and To is the 
time of pericenter passage. The mean motion n = lix /Pb is 
inversely proportional to the time elapsed between two con- 
secutive crossings of the pericenter, i.e. the anomalistic pe- 
riod Pb- In Newtonian mechanics, for two point-like bodies, 
n reduces to the usual Keplerian expression = 2n/P(°>. 
In many binary systems, as in the double pulsar one, the 
period Pb is accurately determined in a phenomenological, 
model-independent way, so that, in principle, it accounts 
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for all the dynamical features of the system, not only those 
coming from the Newtonian point-like terms, within the 
measurement precision. 

The Gauss equation for the variation of the mean 
anomaly, in the case of an entirely radial disturbing accel- 
eration A like eq. (1), is 



dM 2 /r 
—— =n A - 

at na \a 



(1-e 2 ) 



A cos /, 



(4) 



where / is the true anomaly, reckoned from the pericenter. 
Using the eccentric anomaly E, defined as 



M 



esinP, 



(5) 



turns out to be more convenient. The unperturbed Keple- 
rian ellipse, on which the right-hand-side of eq. (4) must be 
evaluated, is 



r = a(l — ecosE) ; 
by using eq. (1) and 
r dM 



(6) 



dE 
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1 — e cos E, 

cos E — e 
1 — e cos E 

eq. (4) becomes 

dE _ n<°> 

dt (1 — e cosP) 



(7) 
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Since Ac 2 /3n 
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from eq. (8) it can be obtained 



(1 — e cosP) 
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which integrated yields that 

P b = P (0 > + P A 

with 

^Ac 2 (7 + 3e 2 ) 



Pa = 



3n(°) c 



(9) 
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2.1. Combining the periastron and the orbital period 

The general relativistic expressions of the post-Kcplerian 
parameters r, s and u> are 



r = T m B , 



is able to test the hypothesis that A ^ 0. From Table 1 it 
turns out 



2- 
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where 



GTtr. 



T Q = 

c° 

and M — tua + ?n-B in units of solar masses. 
By means of 

a = -(xa + xb) 



(13) 



(14) 



and of the equations for r and s it is possible to express 
P' ' and wipn in terms of Pb and of the phenomenologi- 
cally determined Kcplcrian and post-Kcplcrian parameters 
xa,xb, r,s as 



1/2 



3sr 
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[tt(x a + x B )\ {—) s a/ , 
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(15) 



ar A (l - e 2 ) \2n / 

In such a way we can genuinely compare them to Pb and 
Co because they do not contain quantities obtained from the 
third Kepler law and the general rclativistic pcriastron pre- 
cession themselves; moreover, we have expressed the sum 
of the masses entering both p(°> and wipn in terms of r 
and s, thus avoiding any possible reciprocal imprinting be- 
tween the third Kepler law and the periastron rate. At this 
point it is possible to construct 

£• < 16 > 



p 



with 

ACu = Co — wiPNj 

AP = P b -P (0) ; 
note that 

R = R(P h ,XA,x B ,e;u!,r,s). 



(17) 



(18) 



By attributing Ad) and AP to the action of A, not modelled 
into the routines used to fit the PSR J0737-3039A/B timing 
data, it is possible to compare R to 



Pa 



3VT - ^P b r 3 / 2 s 9 / 2 



Pa 4ir 2 {7 + 3e 2 )(x A + x B ) 3 x A /2 



(19) 



and see if eq. (16) and eq. (19) are equal within the errors. 
Note that eq. (19) is independent of A and, by definition, 



Pa = (3.4 ±0.3) x HP 8 s" 



(20) 



Pa is a well determined quantity, different from zero at 
about 11 sigma level. In regard to P we have 



Aw = -0.3 ±2.1 deg yr~\ 



AP = 59 ± 364 s, 



so that 



|P| = (0.3 ±4) x 10 



- 11 s- 2 ; 



(21) 



(22) 



R is compatible with zero in such a way that its range does 
not overlap with the one of Pa: indeed, the upper bound 
on P is three orders of magnitude smaller than the lower 
bound on Pa- Thus, we must conclude that 3 



P^P A . 



(23) 



Concerning the released uncertainties in P and Pa, they 
must be considered as upper bounds since they have been 
conservatively computed by linearly adding the individual 
biased terms due to <5Pb, 5ui, 5e, Sxa, Sx b , Sr, 5s in order to 
account for the existing correlations (Kramer ct al., 2006) 
among them. 

The results of the present study confirm those obtained 
in the Solar System by taking the ratio of the estimated 
corrections to the standard Newtonian/Einsteinian preces- 
sions of the longitude of the perihelia w for different pairs 
of planets (Iorio, 2008). It would be very interesting to de- 
vise analogous tests involving other observablcs (lensing, 
time delay) affected by A as well recently computed in, e.g., 
(Ruggicro, 2007; Sercno, 2008). 

3. The Dvali-Gabadadze-Porrati braneworld 
model 

The approach previously outlined for A can be followed 
also for the DGP braneworld model (Dvali et al., 2000) 
which recently received great attention from an observa- 
tional point of view (Dvali et al., 2003; Iorio, 2008). 

The preliminary confrontations with data so far per- 
formed refer to the perihelia of the Solar System planets. 
Indeed, DGP gravity predicts an extra-precession of the 
pericentre of a test particle (Lue and Starkman, 2003; Io- 
rio, 2005) 



3 / c \ / 13 2 



(24) 



where the signs =F are related to the two different cosmo- 
logical branches of the model and tq is a free-parameter set 



3 In principle, also the 1PN correction to the third Kepler law 
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to about 5 Gpc by Type IA Supernova: data, independent 
of the orbit's scmimajor axis. The predicted precessions of 
about 10 -4 arcsec cy _1 were found to be compatible with 
the estimated corrections to the usual apsidal precessions 
of planets considered one at a time separately (Iorio, 2008) , 
but marginally incompatible with the ratio of them for some 
pairs of inner planets (Iorio, 2007). 

The effects of DGP model on the orbital period is (Iorio, 
2006) 

c_\ a 3 (l-e 2 ) 2 
~ a ) GOT ' 
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"dgp = T^-tt 



(25) 



From eq. (24) and eq. (25) it is possible to construct 



-Rdgp 



^DGP 

-Pdgp 



(26) 



which, expressed in terms of the phcnomenologically deter- 
mined parameters of PSR J0737-3039A/B , becomes 



-Rdgp 



3 (1 _ l§ e 2) Pbr 3/2 s 9/2 



(27) 



227r(l-e 2 )(x A + .T B ) 3 4 /2 ' 

Putting the figures of Table eq. (1) into eq. (27) and com 
puting the uncertainty as done in the case of A yields 



i?DGP = (1.4 ±0.1) x 10~ 7 s~ 2 . 



(28) 



As can be noted, the lower bound of -Rdgp is four orders 
of magnitude larger than the upper bound of R, so that we 
must conclude that, also in this case, 



R 7^ -Rdgp- 



(29) 



The outcome by Iorio (2007) is, thus, confirmed at a much 
more stringent level. 

An analysis of type la supernova (SNe la) data disfavor- 
ing DGP model can be found in (Bcnto et al., 2005). 

4. Conclusions 

In this paper we used the most recent determinations of 
the orbital parameters of the double pulsar binary system 
PSR J0737-3039A/B to perform local tests of two com- 
plementary approaches to the issue of the observed accel- 
eration of the universe: the uniform cosmological constant 
A in the framework of the known general rclativistic laws 
of gravity and the multidimensional braneworld model by 
Dvali, Gabadadze and Porrati which, instead, resorts to a 
modification of the currently known laws of gravity. Since, 
at present, there are no observational evidences for such 
theoretical schemes other than just the cosmological phe- 
nomenon for which they were introduced, it is important 
to put them on the test independently of the cosmological 
acceleration itself. It is worthwhile noting that the results 
for A hold also for any other Hooke-like additional force 
proportional to r. 

To this aim, we considered the phcnomenologically de- 
termined the pcriastron precession uj and the orbital pe- 
riod P h of PSR J0737-3039A/B by contrasting them to the 



predicted 1PN pcriastron rate wipn and the Keplerian pe- 
riod .P^ ) . With such discrepancies we constructed the ra- 
tio R = Aui/AP by finding it compatible with zero: \R\ = 
(0.3±4) x 10~ n s -2 . Then, we compared R to the predicted 
ratios for the effects of A and the DGP gravity-not modeled 
in the pulsar data proccssing-on the pcriastron rate and the 
orbital period by finding R A = (3.4 ± 0.3) x 10~ 8 s~ 2 and 
-Rdgp = (1.4±0.1) x 10~ 7 s~ 2 , respectively. Thus, the out- 
come of such a local test is neatly negative, in agreement 
with other local tests recently performed in the Solar Sys- 
tem by taking the ratio of the non-Newtonian/Einstcinian 
rates of the perihelia for several pairs of planets. 
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Table 1 

Relevant Keplerian and post-Keplcrian parameters of the binary 
system PSR J0737-3039A/B (Kramer et al., 2006). The orbital pe- 
riod P]-, is measured with a precision of 4 X 10~ 6 s. The projected 
semimajor axis is defined as x = (at, c /c) sini, where a^, c is the 
barycentric semimajor axis, c is the speed of light and i is the an- 
gle between the plane of the sky, perpendicular to the line-of-sight, 
and the orbital plane. The eccentricity is e. The best determined 
post-Keplerian parameter is, to date, the periastron rate Co of the 
component A. The phenomenologically determined post-Keplerian 
parameter s, related to the general relativistic Shapiro time delay, 
is equal to sini; we have conservatively quoted the largest error in s 
reported in (Kramer et al., 2006). The other post-Keplerian param- 
eter related to the Shapiro delay, which is used in the text, is r. 

P h (d) x A (s) x B (s) e Co (deg yr _1 ) s r (fjs) 

0.10225156248(5) 1.415032(1) 1.5161(16) 0.0877775(9) 16.89947(68) 0.99974(39) 6.21(33) 



r. 



